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Abstract 

Stokes-Dirac structures are infinite-dimensional Dirac structures defined in terms 
of differential forms on a smooth manifold with boundary. These Dirac structures 
lay down a geometric framework for the formulation of Hamiltonian systems with 
a nonzero boundary energy flow. Simplicial triangulation of the underlaying man- 
ifold leads to the so-called simplicial Dirac structures, discrete analogues of Stokes- 
Dirac structures, and thus provides a natural framework for deriving finite-dimensional 
port-Hamiltonian systems that emulate their infinite-dimensional counterparts. The 
port-Hamiltonian systems defined with respect to Stokes-Dirac and simplicial Dirac 
structures exhibit gauge and a discrete gauge symmetry, respectively. In this paper, 
employing Poisson reduction we offer a unified technique for the symmetry reduction 
of a generalized canonical infinite-dimensional Dirac structure to the Poisson struc- 
ture associated with Stokes-Dirac structures and of a fine-dimensional Dirac structure 
to simplicial Dirac structures. We demonstrate this Poisson scheme on a physical 
example of the vibrating string. 

Keywords: Port-Hamiltonian systems, Poisson structures, Dirac structures, distri- 
buted-parameter systems, symmetry reduction 



1 Introduction 

Geometric structures behind a variety of physical systems stemming from mechanics, elec- 
tromagnetism and chemistry exhibit a remarkable unity enunciated by Dirac structures. 
The open dynamical systems defined with respect to these structures belong to the class 
of so-called port-Hamiltonian systems. These systems arise naturally from the energy- 
based modeling. Apart from offering a geometric content of Hamiltonian systems, Dirac 
structures supply a framework for modeling port-Hamiltonian systems as interconnected 
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and constrained systems. From a network-modeling perspective, this means that port- 
Hamiltonian systems can be reticulated into a set of energy-storing elements, a set of 
energy-dissipating elements, and a set of energy port by which the interconnection of 
these blocks and environment is modeled. It is well-known that such a modeling strategy 
also utilizes control synthesis for these systems. 

The port-Hamiltonian formalism transcends the lumped-parameter scenario and has 
been successfully applied to study of a number of distributed-parameter systems [7, 4] . The 
centrepiece of the efforts concerning infinite-dimensional case is the Stokes-Dirac structure. 
The canonical Stokes-Dirac structure is an infinite-dimensional Dirac structure defined 
in terms of differential forms on a smooth manifold with boundary. The Hamiltonian 
equations associated to this Dirac structure allow for non-zero energy exchange through 
the boundary. 

Although the differential operator in the Stokes-Dirac structure, in the presence of 
nonzero boundary conditions, is not skew-symetric, it is possible to associate a (pseudo- 
)Poisson structure to the Stokes-Dirac structure [7]. In the absence of algebraic constraints, 
the Stokes-Dirac structure specializes to a Poisson structure [3], and as such it can be 
derived through symmetry reduction from a canonical Dirac structure on the phase space 
[10]. How to conduct this reduction for the Poisson structure associated to the Stokes- 
Dirac structure on a manifold with boundary is the central theme of this paper. 

Contribution and outline. This paper very closely follows [10]. The reduction 
scheme we are dealing with is the one from [10], the only difference being in that we con- 
sider slightly augmented spaces in order to account for the behaviours associated with the 
boundary. The perspective as well as the notation in Section 2 are taken verbatim from 
[10], but now for the generalized Dirac structures that allow for the formulation of open 
Hamiltonian systems. The proposed Poisson reduction is firstly applied in the reduction 
of a generalized cannonical Dirac structure to the Poisson structure associated with the 
Stokes-Dirac structure. In the context of dynamics, the canonical port-Hamiltonian sys- 
tems are those defined as in [5, 6], now only in the context of differential forms, while the 
reduced port-Hamiltonian systems are exactly those presented in [7] . In the final section we 
demonstrate how this reduction applies to the Poisson reduction of the port-Hamiltonian 
systems on discrete manifolds [8, 9]. 

2 Dirac structures and reduction 

Dirac structures were originally developed in [2, 3] as a generalization of symplectic and 
Poisson structures. The formalism of Dirac structure was employed as the geometric 
notion underpinning generalized power-conserving interconnections and thus allowing the 
Hamiltonian formulation of interconnected and constrained dynamical systems. 
Let Q be a manifold and define a pairing on TQ ®T*Q given by 



For a subspace D of TQ © T*Q, we define the orthogonal complement D as the space 
of all (v,a) such that (((v, a), (w, /?))) = for all (w,/3). A Dirac structure is then a 
subbundle D of TQ ®T*Q which satisfies D = D^. 
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The notion of Dirac structures just entertained is suitable for the formulation of closed 
Hamiltonian systems, however, our aim is a treatment of open Hamiltonian systems in 
such a way that some of the external variables remain free port variables. 

Let F be a linear vector space of external flows, with dual the space F* of external 
efforts. We deal with Dirac structures on the product space Q x F. The pairing on 
(TQ x F) © (T*Q x F) is given by 

((v,f),(a,e)),((wj~),((3,e)))) 

1 / \ (2- 1 ) 

= -{a(w) + e(f)+f3(v) + e(f)). 

A generalized Dirac structure D is a subbundle of (TQ x F) © (T*Q x F) which is 
maximally isotropic under (2.1). 

Canonical Dirac structure on TQ © T*Q is considered to be a symplectic structure. 
However, in this paper we shall deal with slightly different canonical Dirac structures. To 
that end, let the map (J : T*Q x F* — > TQ x F induces the Poisson structure on TQ x F. 
The graph of fl given by 

Dt*QxF* ■= {(tt(a, e), (a, e)) : a G T*Q , e G F*} (2.2) 

is a Dirac structure. If the mapping fl is symplectic on TQ, that is if j)(a, 0) = implies 
a = 0, the Dirac structure (2.2) is the generalized canonical Dirac structure. 

There is a number of techniques for symmetry reduction of Dirac structures [1, 11]. The 
reduction considered in this paper is the Poisson reduction from [10]. For that purpose, let 
G be a Lie group which acts on Q from the right and assume that the quotient space Q /G 
is again a manifold. Denote the action of g G G on q G Q by q-g and the induced actions of 
g G G on TQ x F and T*Q x F* by (v, f)-g and (a, e)-g, for u G TQ, / G F, a G T*Q, and 
e G F*. The action on the T*Q x F* is defined by {(a, e) ■ g, (v, /)) = ((a, e), (u , /) • g^ 1 )- 
In what follows, we will focus mostly on the reduced cotangent bundle (T*Q x F*)/G. In 
this paper, we will deal with the space denoted by T*Q/G x F*. 

Consider now the canonical Dirac structure on T*QxF*. Let (J : T*QxF* — > TQxFbe 
the map (2.2) used in the definition of Dt*qxF* ■ The reduced Dirac structure -Dt*q/GxF* 
on T*Q/G x F* can now be described as the graph of a reduced map [ft] : T*(T*Q/G x 
F*) -> T(T*Q/G x F*) defined as follows. 

Let 7rc : T*Q x F* — > T*Q/G x F* be the quotient map and consider an element 
(p,ir,Pb) in T*Q x F*. The tangent map of ttg at (/?, 7r,p&) is denoted by T^ p7Tpb ^-KG ■ 
T^ Pb) {T*QxF*) -4 T(p,n,p b ) (T*Q/G x F*) , and its dual by Z^tto : Ko^fTQ/Gx 
F*) — > T*, *.(T*Q x F*). The reduced map [jj] now fits into the following extended 
commutative diagram 



r { ;.,p 6) ( T ^ >< n — - — -?w 6) (t*q x F *) ( 2 - 3 ) 
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3 Constant Stokes-Dirac structures 



Throughout this paper, let M be an oriented n-dimensional smooth manifold with a 
smooth (n — l)-dimensional boundary dM endowed with the induced orientation, repre- 
senting the space of spatial variables. By Q, k (M), k = 0,1,..., n, denote the space of 
exterior /c-forms on M, and by Q k (dM), k = 0, 1, . . . , n — 1, the space of fe-forms on dM. 
A natural non-degenerative pairing between p G f2 fc (M) and a G Q n ~ k (M) is given by 
(a\p) = f M o~ A p. Likewise, the pairing on the boundary dM between p G Q k (dM) and 
a € O re - fe - 1 (5M) is given by (a\p) = f gM a Np [7]. 

3.1 Stokes-Dirac structure 

For any pair p, (/ of positive integers satisfying p + q = n + 1, define the flow and effort 
linear spaces by 

3 M =Q p {M) x fl q (M) x O n ~ p (dM) 

£ M =n n - p (M) x O n -«(M) x ST^dM) . 

The bilinear form on the product space 3p i(? x £ Pj(? is 

el ^1 ^1 „1 „1 „1\ f 2 2 2 2 



(({fpifq^fb^pi^q^b)iifpifqifbi^pi^q^b"))) 

= / elAfZ + elAfZ + elAf' + elAft (3-1) 
Jm 

+ [ elAft + elAfl. 
JdM 



Theorem 3.1 (Stokes-Dirac structure [7]). Given linear spaces Jp^ q and £ p ,q, and the 
bilinear form ((,)), define the following linear subspace T> of'Jp^ x E pq 



2) — { (/pj /qj fb-i %>) ^-q-, &b) € 3~p,q X £p j( 

/ P U f (-ir+M \ / e p 

/J V d J\e g J ' (3.2) 

A A = ( tr \ / e p 

e b J \ -(-iy-«tr J\e q 

where d is the exterior derivative and tr stands for a trace on the boundary dM. Then 
T> = D , that is, T> is a Dirac structure. 

It is possible to associate a Poisson structure to the Stokes-Dirac structure D. Here 
we just sketch the essence and refer the reader to [7]. 

The space of admissible efforts is £ a dm := {e G £ P) g|3/ G 3~ Pi(3 such that (/, e) G 
D}. The set of admissible mappings %sd m := {& : 3~ Pi(? -4 M|Va G 3~ P , q , 3e(k,a) G 
£adm such that for Va G 3p ig k(a + 9a) = /c(a) + (e(k,a)\da) + 0(<9a)}. The set 3C a dm 
consists of those functions k : fF(M) x il 9 (M) x O n_p (dM) — )• M whose derivatives 
fffe(^) = (5pk(z),5 q k(z),S b k(z)) G O n -P(M) x O re -«(M) x n n ~i{dM) satisfy = 
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-(— l) n q tr(5gk(z)). The Poisson bracket on 3C a d m is given as 

{k\ k 2 }^ = [ ((dpk 1 ) A (-l) r d((5 q k 2 ) + (Sqk 1 ) A d(6 p k 2 )) 

J M 

((-1)^(5^) A (5 p k 2 )). 



8M 

Using Stokes' theorem, it follows that the bracket is skew-symmetric and that it satisfies 
the Jacobi identity: {{k 1 , k 2 }?,, k 3 }v + {{^2, k 3 }®^ 1 }® + {{k 3 , k 1 }®, k 2 }?> = for all 

k (E 3<Cadm' 

In this paper we will exclusively be dealing with Poisson and associated Poisson struc- 
tures. 

3.2 Simplicial Dirac structures 

In the discrete setting, the smooth manifold M is replaced by an n-dimensional well- 
centered oriented manifold- like simplicial complex K [8, 9]. The flow and the effort spaces 
will be the spaces of complementary primal and dual forms. The elements of these two 
spaces are paired via the discrete primal-dual wedge product. Let 

9* = ngfaiir) x n d (iq x n n d - p (d(K)) 



p,i 



The primal-dual wedge product ensures a bijective relation between the primal and 
dual forms, between the flows and efforts. A natural discrete mirror of the bilinear form 
(3.1) is a symmetric pairing on the product space S~ q x £^ defined by 

((( fpi fqi fb i e pi &\i &b)i (fpi fqi fbi e pi ^qi^b)))d 

= <ej A f 2 + e]A f 2 + e 2 p A f l p + e\ A f q , K) 
+ (elAf 2 + e 2 Af b \dK). 

A discrete analogue of the Stokes-Dirac structure is the finite-dimensional Dirac structure 
constructed in the following theorem [8]. 

Theorem 3.2 (Simplicial Dirac structure [8]). Given linear spaces 3"p q and &p q , and the 
bilinear form ((,))d- The linear subspace D d C x E pq defined by 

= {{fpi fqi fbi e pj &qi (Zb) & ^p^q x ^p,q\ 

!,\_( (-irdTV* V(-irf d r^ 6 , (3-4) 



f q J V dn_P 
fb = (-l) p tr""%}, 

with r = pq + 1, is a Dirac structure with respect to the pairing ((, )) d . 
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The operators d n p is the discrete exterior operator mapping $7^ P (K) to tt q d (K), and 
d™~ q is the dual discrete exterior derivative. Note that since d^~ q = (— l) q (d n ~ p ) T and 
d^~ q = (— l) n ~ p (tr n ~ p ) T , the structure (3.4) is in fact a Poisson structure on the state 
space n p d (*iK) x ^(if). 

The simplicial Dirac structure (3.4) is used as terminus a quo for the geometric for- 
mulation of spatially discrete port-Hamiltonian systems [9] . 



4 Reduction of Stokes-Dirac structure 

The configuration manifold is a vector space Q := tt k (M) with the tangent bundle TQ = 
Q x Q and the cotangent bundle T*Q = Q X Q*, where Q* = Q n ~ k (M). The space of 
the boundary flows F will be an admissible subset of VL n ~ k ~ 1 (dM), while the space of the 
boundary efforts is E := F* = Q k (dM). 

The tangent bundle T(T*Q x F*) is isomorphic to (Q x Q* x F*) x (Q x Q* x F*), 
with a typical element denoted by (p, tt, p b , p, n, pb), while T*(T*Q x F*) = (Q x Q* x 
F*) x (Q* x Q x F), with a typical element denoted by (p, tt, pb, e p , e n , e b ). For the duality 
pairing between T(T*Q x F*) and T*(T*Q x F*) we chose 

((p, it, p b , e p , e n , e b ), (p, tt, p b ,p, tt, p b )) 

f f (4- 1 ) 

= / (e p A p + e n A 77 ) + / (e& A p b + e b A tr p) . 

Jm JdM 
The choice for this non-degenerate pairing will become clear later on. 



4.1 The symmetry group 

Let G be an Abelian group of (k — l)-forms. For any a G G and p G Q, the group G 
action on Q is 

p-a = p + da. (4.2) 

This action of gauge group lifts to TQ x F and T*Q x F* as (p, p, e&) • a = (p + da, p, e&) 
and (p, 77, p fe ) - a = (p + da,77, p b ) for a <G G, (p,p,e 6 ) € TQ x F and (p,ir,p b ) G T*Q x F*. 

The elements of Q/G are equivalence classes [p] of A;- forms up to exact forms, so 
that the exterior differential determines a well-defined map from Q/G to df2 fc , given by 
[p] i— )> dp. If the k-th. cohomology of M vanishes, we have Q/G = dQ k . Consequently, the 
quotient (T*Q/G x F*) is isomorphic to Q/G x Q* x F*, or explicitly 

(T*Q x F*)/G = dn k (M) x n n ~ k (M) x n k (dM). 

The quotient map denoted as ttq ■ T*Q x F* — > (T*Q)/G x F* is given by 

?Tg(P, 77, p b ) = (dp, 77, p b ). (4.3) 

Let a representative element of T*Q/G x F* be (p, 7f,p&), with p G d$7 fc (M), 77 G 
O n " fc (M) andp b G Q k {dM). Elements of T(T*Q/GxF*) will be denoted by (p, tt, p 6 , p", #, p 6 ), 
while the elements of T*(T*Q/G x F*) will be denoted by (p, 77, p&, e p , e^, e&). For the du- 
ality pairing, we use 

{(p,Tt,p b ,e p ,e n ,e~b), (p, 77, p b , p,Tf,pb)) = \ (e p A p + A 77) + / e b A p b . (4.4) 
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Whenever the base point (p, n, pb) is clear from the context, we will denote (p, ir, pb, p, tt, pb) 
simply by (p\ 

Pb): and similarly for (p, tt, pb, e p , e n , e b ). 
4.2 The reduced Dirac structure 

The generalized canonical Dirac structure is a Poisson structure induced by the linear 
maping jj : T*(T*Q x F*) -> T(T*Q x F*) given by 

$(p,n,p b ,e p ,e n ,e b ) = (p, tt, p b , e n , -(-l) fc(n_fc) e p , -tr e n ). (4.5) 

In order to obtain the reduced Poisson structure from the canonical Dirac structure 
(4.5), we need to specify what are the operators Ttt g and T*tt g in the diagram (2.3). The 
space F is the set of admissible forms Q n ~ k ~ l (dM) that are the traces of (dO fc )*, as will 
be made clear in Lemma 4.1. Consider an element (p,ir,p b ) € T*Q x F*, and we recall 
that TT G (p,TT,p b ) = (dp,TT,p b ). Let T (Pi7rjPb) TT G : T^^ Pb) (T*Q x F*) -4 T {dp ^ )Pb) (T*Q/G x 
F*) be the tangent map to tt g at (p,n,p b ) and consider the adjoint map T* p ^ p A^G '■ 
T*^, Pb) (T*QlG x F*) -4 T ( ;^ jP6) (T*Q x F*). 

Lemma 4.1. T/ie tangent and cotangent maps T^p^^itQ and T^^^sWg are given by 

T(p,-K,p b )KG{P, 7T, Pfe, p, 7T, p h ) = (dp, 7T, p b , dp, TT, p b ) (4.6) 



and 



T (p,n, Pb ) 7r G{dp,7r,p b ,e p ,e 7r ,-{-l) n tre p ) 
= (p, tt, p b , (-l) n - k de p , e n , -(-lf^tr e p ). 



Proof. The expression (4.6) for Ti p ^ pb ^TTG follows from (6.3). To prove (4.7), we let 
(p,n,p b ) G T (/VriPb) (T*Q x F*) and consider 



T^n^Giep, e T , -(-l) n_fc e p ), (p, tt, p 6 ) 
= ((e p , e^, -(-l) n ~ fc tr e p ), T (/VrjP()) Tr G (p, tt, p b )^> 
= ((e P ,e 7r ,-(-l) n - fc tre p ),(dp,7r,p 6 )^) . 

Applying Stokes' theorem, we have 

{{e pi e n , -(-l) n - k tre p ), (dp, tt, p 6 )^ 

= / (e p A dp + e„- A tt) 

+ / (-(-l)"- fc tre p Ap b -(-ir- fc tre p Atrp) 
Jm 



((-l) n_fc de p A p + A tt) - / (-l) n - fc tr e p A tr p . 
Thus, r ( ;^ pb) 7T G (e p ,^,-(-l)"- fc tre p ) = ((-l)»-*de„ e*, -(-l^tr e p ). □ 
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As in the case of a boundaryless manifold [10], the reduced Poisson structure in (2.3) 
is given by 

M(dp,n,p b ) = T { p,TT,p b )KG ° It O T^^TTG 

for all (dp,7r,p 6 ) G T*Q/G x F*. 

Theorem 4.2. T/ie reduced Poisson structure is given by 

[tJ](e p ,e^,-(-l) n - fc tre p ) = (de*, _(_l)»(*+i) dgp) -tr e*). (4.8) 

Relation to the Stokes-Dirac structure. The matrix form of the reduced Poisson 
structure is 

d 0\ / e p \ 

_ 1 )n(fc+i) d o o J e w . (4.9) 

-tr 0/ \(-l) n -*tre p / 

The sign convention in (4.9) and [7] is not the same. To match the signs we introduce new 
flow variables f p , f q , fb and effort variables e p , e q , eb defined as e p = e p , e q = (—l) r e n , f p = 
Pi fq = ( — l) Tl ( fc + 1 ) +1 7f , f b = —(—l) r p b: where p = k + 1, q = n — k, and r = pq + 1. With 
this choice of signs, (4.9) becomes 











w 





7A = (0 (-l) r d\ (e p 
JJ \d 
fb = tr e q 
{-l) n - k tre p = e b . 



(4.10) 



Here, it is important to point out that the boundary effort eb, unlike in the case of the 
Stokes-Dirac structure, does not follow from the associate Poisson structure, but rather 
belongs to the set of admissible derivatives of the flow restricted to the boundary. 



5 Symmetry in port-Hamiltonian systems 

Let t i—)- {pip, oijc) £ Q k (M) x Q n ~ k (M) be a time function, and let the Hamiltonian be 

H(a p ,a jr )= / JC(da p ,an). 
Jm 

It follows that at any time instance i G R 

dH _ r 5H da p 5H ^ da^ f d'K da p 
dt J M 5a p dt Sa^ dt J gM d(da p ) dt 

The differential forms , represent the generalized velocities of the energy variables 
ctp,a.-fr. The connection with the canonical Dirac structure is made by setting the flows 

da p da-jr 
P = ~~dt 7F= 9T' 



S 



and the efforts 

SH 5H 

oa p oa^ 

The canonical distributed- parameter port-Hamilto-nian system on an ?i-dim- 
ensional manifofd, with the state space £l k (M) x Q n ~ k (M), the Hamiltonian H and the 
canonical Dirac structure (4.5), is given as 

(ii) - u t) (t ) ■ 

Proposition 5.1. For the port- Hamiltonian system (5.1) the following property 

-TT = / e 6 A / 6 
at JaAf 

expresses the fact that the increase in energy on the domain M is equal to the power 
supplied to the system through the boundary dM . 



5.1 The reduced port-Hamiltonian systems 

The Hamiltonian H is invariant if a spatially independent /c-form is added to «p, thus 
the Poisson reduction is applicable. Let the reduced field be <5p := dap, then the reduced 
Hamiltonian is 

H(a p ,ait)= \ ^(a^aj,). 
J M 

The port-Hamiltonian system with respect to the reduced Poisson structure is 




d 



f MV 
I 5a e 



-(-l)"Md o) [jl 


-(-l) n " fc tr 




(5.2) 



This is precisely the port-Hamiltonian system given in [7]. 

We will show how the general considerations of the reduction of port-Hamiltonian 
systems apply to a physical example of the vibrating string. 



5.2 Vibrating string 

Consider an elastic string of length I, elasticity modulus T, and mass density /i, subject to 
traction forces at its ends. The underlying manifold is the segment M = [0, 1] C M, with 
coordinate z. 

Under the assumption of linear elasticity, the Hamiltonian is given by 

H(u,p)= I 3t(u,p) = - [ {^r 1 pA*p + TduA*du), 
Jm 2 J M 
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where p 6 Q l (M) is the momentum conjugate to the displacement u G f2°(M), and * is 
the Hodge star. 

The canonical Hamiltonian equations are 

/ 8u\ / q ^\ / SH 



w 



fb = tr(* M - 1 p) (5.3) 



e b = tr 



U(d«)J 



or component- wise 



du _i 
ft = *"~ P 

^ = d(*Tdn) 

/& = tr(*/^~ 1 p) 
ej, = tr(*Tdu) . 

The Hamiltonian formulation (5.3) is identical to the formulation of the heavy chain system 
in [6]. 

The energy balance for the vibrating string is 

dH _ f 5H du 5H dp f d'K du 
dt J M 5u dt 5p dt Jqm d(du) dt 



■ / -d(*Tdu) A*u 1 p + */z 1 pAd(*Tdu) 
JM 

+ / *u~ x p A *Tdu 
JdM 

■ / *u~ l p A*T du = / e b Af b . 
JdM JdM 



The Hamiltonian is invariant if a time function is added to u. The potential energy 
can be expressed in terms of the strain a = du so that the reduced Hamiltonian is given 
by 

H r (a,p) = / %r(u,p) = - / A *p + Ta A *a) . 

Jm 2 J M 

The Hamiltonian equations of the vibrating string now read as 



fb) = (0 tr 
eJ Ur 



SH r 
Sa 

SH r 
6p 



(5.4) 



These are the equations that correspond to the formulation of the vibration string system 
with respect to the Stokes-Dirac structure as is given in [7]. 
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6 Symmetry reduction in discrete setting 

In the discrete world, the configuration space is the set of primal discrete forms Q = Q k (K) 
with the dual Q* = n n - k (*iK). The space of the boundary efforts is E = F* = n k (d(K)), 
and the space of the boundary flows is F = Q, n ~ k ~ 1 (d(*K)). 

For the duality pairing between T(T*Q x F*) and T*(T*Q x F*) we choose 



{(p,Tr,Pb,e p ,e nj e b ), (p, vr, p b , p, vr, p' fe )) = / (e p A p + e n A vr) + / e b A p b , (6.1) 

JM JdM 

where A is the primal-dual wedge product. 

The generalized canonical Dirac structure is a Poisson structure induced by the linear 
mapping (j : T*(T*Q x F*) T(T*Q x F*) given by 

Jj(p, 7T, p 6 , e p , e*, e 6 ) = (p, vr, p b , e n , -{-l) k{n ~ k) {e p + dj* - * -1 ^), -tr k e n ). (6.2) 
The group G that acts on Q is described by the following action 

lk— 1 

a • p = p + a a 

for a € G and p £ Q, where d fc_1 is the discrete exterior derivative. 

The quotient is (T*Q/G x F*) = d k il k {K) x O n " fc (* iJ FC) x VL k (d{K)). 

As in the continuous setting, the quotient map denoted as txq : T*QxF* —¥ (T*Q) /Gx 
F* is given by 

ttg(p, tt, Pb) = (d fc p, vr, p fe ). (6.3) 
For the duality pairing between T*(T*Q/G x F*) and T(T*Q/G x F*), we take 



((p,Tr,p b ,ep,e n ,e b ),(p,n,p b ,p,Tr,p b )) = (e p A p + e n A tt) + / e b A p b . 

JM JdM 

As before, whenever the base point (p, vf, p~b) is clear, we will denote (p, vr, p~ b , p, f, p~b) simply 
by (p,^r,p b ), and similarly for (p, 7f, pi, e p , B n , e b ). 

Lemma 6.1. The tangent and cotangent maps T(p^ jPb )irG and ^( p7rpi) ) 7r G are given by 

T(p,n,p b )KG{P, TT, Pfo, P, TT , p b ) = (d fc p, TT, p b , d fc p, TT, p b ) (6.4) 

and 

T ( P ,TT, Pb ) 7T G(d k p, vr, pb, e p , e^, e b ) = (p, vr, p 6 , (-l) n " fc d™~ fc ~ 1 e p , e^, e 6 ). (6.5) 

Theorem 6.2 (Reduced simplicial Dirac structure). The reduced simplicial Poisson struc- 
ture is given by [«](e p , e v , e b ) = (d fc e w , -(-l^+^-lf-^r^+dr* -1 ^), -tr*^). 
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Port-Hamiltonian systems on a simplicial complex 

The canonical port-Hamiltonian system with respect to the canonical Dirac structure is 

da p _ dH 

(-!)*<"-*> ( ^(a^a^ + d^e,) (6.6) 



dt \da p j 

jL dH . . 

p b = -tr ^— (Q!^,Q4r) 

The rank of the underlying Poison structure is the rank of the symplectic phase space 
n k (K) x Q, n - k (*iK). 

The canonical Hamiltonian (ap,a^) i— > H(a p ,a^) can be expressed as 

H{a p ,ajt) := i? r (d fc a p ,a^-) . (6.7) 

The reduced port-Hamiltonian equations assume the following form 

= =d^(« p ,a.) = d^K,a.) 

ot ot oa-jr da^ 



p b = -tr ^— (ap,a,v) = -tr ^— - (a^, a*) 



This is precisely the port-Hamiltonian system on a simplicial manifold as presented in 
[8, 9]. 



7 Final remark 

This paper addresses the issue of the symmetry reduction of the generalized canonical 
Dirac structure to the Poisson structure associated with the Stokes-Dirac structure. The 
open avenue for the future work is to find a reduction procedure that would directly lead 
to the Stokes-Dirac structure. 
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